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Abstract 

A simple formula for the scattering of wave packets from a square well at long times 
is derived. The expression shows that the phenomenon of wave packet diffraction in 
space and time exists in three dimensions also. An experiment for the verification of 
the phenomenon is described and simulated numerically. 
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1 Introduction 

In a previous work|l|, it was shown analytically that there exist a broad class of phenomena 
occurring in nonrelativistic and relativistic quantum wave packet potential scattering that 
behave as time dependent persistent diffraction patterns. The patterns are produced by a 
time independent well or barrier. Differing form the diffraction in time[0] and diffraction in 
space and time|§ phenomena, that apply to the time dependent opening of slits with plane 
monochromatic waves, the pattern does not decay exponentially with time. 

Diffraction in time [[J], and diffraction in space and time|| are very valuable tools to inves- 
tigate time dependent predictions of quantum mechanics. Recent measurements of atomic 
wave diffraction support the theoretical prediction of the diffraction in time process. The 
phenomenon of wave packet diffraction in space and time is still untested. We will propose 
below an experiment for its verification together with a numerical simulation of the expected 
results. 

In ref.|TJ a closed analytical formula for the wave function of a packet scattering in 
one dimension off a time independent potential was derived. This work was preceded by 
numerical investigations of the one and two-dimensional cases. [f|, || It was found analytically 
and numerically that, a diffractive structure that travels in time, exists for all packets, but, 
it survives only for packets that are initially narrower than a value related to the potential 
effective extension. For wider packets, the peak structure merges into a single peak. 

In the present work we develop an analytical formula for the three dimensional case and 
show that the same behavior is present. In section 3 we also treat the one dimensional Dirac 
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case in an analytical manner. Section 4 proposes an experiment and simulates its expected 
results. 

2 Wave packet diffraction in three dimensions 

Consider a gaussian wave packet, 



0(r, t) = e vo 

. ^ / — * — * \ (f-r ) 2 
y = i q • (r - r ) ^ 2 (1) 

with q , the initial average momentum of the packet, r , the initial location and a the 
packet width, impinging on a spherically symmetric square well 



V(x) = -V Q{w - r) (2) 

where, w is the well width, and Vb, the depth, and 9, the step function. The results are 
not dependent on the choice of wave packet [Q], however, the Gaussian profile facilitates the 
integrations. The total wave function with outgoing boundary condition reads 



ip + (r,t) = ip in (r , t) + ip scatt (r , t) 

= J d 3 k /(k, q ) (^n(k, r, t) + ?pscatt(K r, t)) (3) 



where 
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^n(k,?,t) = e lk ' ? 

</w(k,r,t) = ^£(2L + 1) «n(<fr) P £ (*-k) 
ft r L 

/(k, q ) = e~£-^ 2 C *<>-* & * (4) 
The incoming wave may be found explicitly to be 



^ n Gt/(2m) + a 2 ) 

(r - r ) 2 - 4 i (r - r ) ■ g a 2 - 4<x 4 gg 
yi " zt/(2m) + a 2 (5j 

Performing the angular integral over k, the scattered wave becomes 



ip 'scatt(r, t) = J k 2 dk <5>(k,t) 

$(ft, t) = ^—^ Y&L + 1) e iSL stn(5 L ) I L e~ a ' ( fc2 +«oH 
k r ^ 

(6) 

The expression for is quite involved. However, for L = it simplifies considerably. 
Fortunately, this is the only contribution we need at long times (except for very specific 
cases for which the phase shift of higher partial waves dominates, see below) . Fot t — > oo the 
integral above oscillates wildly except for values of k that give an extremal phase |J. It is 
easy to show that in this case, only extremely small values of k enter the integral. However, 
for very small values of k, the phase shift reads 
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(fc W ) 2 ^ g - / 

tan{6l) = -(2 /-!)!! (2/ + l)!!^TTTT 



where x = k' w = yk 2 + 2 m \Vq\ w, and ji, denotes the spherical Bessel function of 
order ]. 

It is then clear, that, except for very pathological cases only the lowest partial wave 

w x l 



matters. Higher partial waves are suppressed by factors of the form , . For 

\yjit/(2m) +a 2 J 

very long times we can then approximate the scattered wave by taking only the L = term 
in eq.(||). For this case II becomes 



I = An j (k d) 

d = |r + 2i(T 2 qo| (8) 

where the absolute value pertains to the vectorial character of the variable inside the 
vertical bars, and not to the complex number. Therefore, jo is really a complicated mixture 
of Bessel and Hankel functions. 

Inserting eqs. (|8],|7|) i n eq. (0) and recalling the rules for the integration of the error function 
for complex variables, we find the transparent result 



^scatt 
2/2 



7T 



it/ (2m) +a 2 
(r + df 



4 {it /(2m) 



r + d 



+ A 



(9) 
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where ln(X) = — — -, with z\ defined in eq.(0). 

a (z + 1) 

The incoming and scattered wave interfere in the total wave. Noting that the scattered 
wave depends on the absolute value of r, while the incoming wave depends on r itself, it is 
clear that there will be a totally different behavior at forward and backward angles. The 
backwards scattered full wave for long distances and times becomes 

/ jy \ 2 777, 7* 

* = "*V - ^ + i± r 2 + A/2 (10) 

This is a diffraction pattern that travels in time as found for the one dimensional casejl]]. 
Again the pattern gets blurred and forms a single peak when the imaginary part of the 
argument of the sin function is large. As explained in ref. ||], this amounts to the condition 
(r » For narrow enough packets the pattern persists to infinity. In the forward 

direction there is no interference as may be easily found by examining eq.(|9|). Instead of a 
sin function, we find there a pure exponential. Only tiny interferences survive due to small 
terms. In some sense the process should be more properly called, interference in space and 
time, because it is due to the superposition of incoming spread packet and the scattered 
wave. 

Numerical calculations of the two dimensional case|| indeed support these findings. 
Moreover, the results are independent of the impact parameter provided the initial packet 
is located at a much longer distance along the direction of q as compared to the transverse 
direction, as found numerically in ref.||. Also the pattern does not depend strongly on the 
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initial momentum. The formula of eq.(|T(]) may be used to find the angular distribution, 
the scattering and total cross sections, etc. Such expressions do not seem to exist in the 
literature for packets, even for the simplest square well case dealt with here. 

As mentioned above the above treatment took advantage of the properties of Gaussian 
packets and the simplicity of the square well. However, they are general. One has only 
to replace the Gaussian form factor by the appropriate one in case and the phase shift for 
the square well by the one corresponding to the specific choice of potential. The case of 
long range potentials deserves however some care. In refs.fl], [5], |6| the connection to the 
excitation of a half-bound state, inside the well, Levinson's theorem and the phenomena 
of diffraction in time and diffraction in space and time were discussed. The same applies 
without modification to the present case. 

It is also evident, that the sign of the potential is irrelevant. A barrier may do too.(l| 
This is in line with optical diffraction of slits and around objects. 

3 Wave packet scattering for the Dirac equation 

Consider the one-dimensional Dirac equation with a scalar potential S(x) as given by eq.(^) 

^ = (11) 
For the initial packet we use a minimal uncertainty relativist ically invariant wave packet 



«7 M dp + (m + S(x)) 
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u 

V 



= J dk e l k ( x ~ x »)~ 2 °" 2 ( E E o~k qo-m 2 ) 

-I 



E + m 



J k(x—xo)—2 a 2 (E Eo—k go—m 2 ) 



(12) 



with E = \/k 2 + m 2 , E = \J q$ + m 2 , and we have taken 70 = <r z , 7i = &z Cx for the 
Dirac matrices, where a denotes the corresponding Pauli matrix. 

The stationary wave function inside and outside the well for fixed momentum is given by 

4>(k, x, t) = 



(f>i{x < —w, k, t) = 

4> 2 {x < —w, k, t) = 

<pi(—w<x<w,k,t) = 

4>2(—u;<x<w,k,t) = 

4>\{w < x, k, t) = 

2 (^ < x, k, i) = 



(e ikx + B e- { k x ) 
i k 



(e l k x - B e~ l k x ) 



E + m 

(C e* k ' x + D e _i k ' x ) 

i k' 
E + m* 

p £ i k x 

i k 



(C e ik ' x - D e~ l K x ) 



i k' x\ 



E + m 



F e 



i k x 



(13) 



with 



B = 



(l-9 2 M-el) 
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c 

D 
F 



-2 ei e 4 jl + g) 
A 

2 ei e 3 (1 - #) 



A 



-Ag_ 
A 



(14) 



where ei = e* fc ™, e 2 = e~ l kw , e 3 = e* k ' w , e 4 = e"* k ' w , g = 7 ^= ' ^ , A 



fc' (E + m) 



(E + m*)' 



e 2 (1 - ^) 2 - ej (1 + (?) 2 , fc' = \J E 2 - (m*) 2 , m* = m - H-Q 
The solution to the scattering of the packet becomes *f?(x,t) 



( \ 

U(x,t) 

v V(x,t) t 



U(x,t) 
V{x,t) 



dk Mk,x,t) e~ l k Xo ~ 2 a2 {E Eo ~ k -?o-™ 2 )-* e t 
dk (t> 2 (k,x,t) e~ l k xo ~ 2 ° 2 {E Eo ~ k w-™ 2 )" 1 E 1 



(15) 



As in section |2] we consider the long time behavior of the wave. Here too, the reflected 
wave will receive contributions mainly from low values of k. Inspection of the reflection 
parameter B in eq.([T3|), we find that for k —>■ 0, B — > — 1. The same happened in the 
nonrelativistic treatment of ffl. For low initial momenta, go << rn,or v << 1, the expo- 
nential of eq.flT5|) becomes « e~ ik m {q °- k)2 . So except for the factor e l m *, 
that is independent of k, we are back at the nonrelativistic case for the upper component, 
and at the same time the lower component vanishes. Hence, the Dirac case in one dimension 
is exactly analogous to the nonrelativistic case for low average momenta of the initial wave 
packet. 



1 For a vector potential the * is in the energy instead of the mass 
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In the ultrarelativistic limit m — > 0, the upper and lower components are identical up to 
a factor of i. The incoming and reflected waves do not spread. Hence, there is no interference 
between both. If we first find the wave in the backward direction as in ref . QXj , and then take 
the limit of a small mass, we find that the pattern disappears, as the arguments of the sin 
function tend to zero. For t — > oo, the backwards scattered wave becomes a receding wave 
packet located around x = —xq —t — w, without any diffraction. The incoming packet is no 
longer present in the region behind the well. 

In summary, the diffraction in space and time with wave packets exists in the relativis- 
tic regime, but gradually diminishes as the mass decreases, until is eventually washed out 
completely for the massless case, as was found for the diffraction in time process with the 
Klein-Gordon equation. f2[ 

4 Numerical simulation of a suggested experiment 

Present advances in bose traps permit the easy handling of atoms at very low temperature 
and velocities. At such energies the atoms move as if they were a wave packet provided the 
random agitation due to thermal effects is not as crucial. However, still the large amount of 
interactions between the particles and of them with the environment causes a decoherence 
of the wave in quite a short time. 

Instead of resorting to a cold bose gas, we will then focus on a drop of liquid Helium. 
Although the experimental details are beyond the expertise of the present theoretical work 
we will describe a hypothetical setup in the hope that it is not too far away from reality. 
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Consider a drop of liquid helium, of around 1 cm 3 in volume lying at a corner of Cesium 
coated plates perpendicular to each other. Cesium is needed in order to prevent wettingQ. 
The helium, plate and environment should be at a temperature below the A point, although 
it would be interesting to see the effect of normal to superffuid component of the Helium on 
the wave packet behavior, as a tool to study decoherence. Put a large area detector behind 
the drop behind the plate. The number of particles measured as a function of time at a fixed 
position would then be determined by the absolute value of the wave function at that point. 
The particles are repelled by the wall and the packet spreads at the same time. Figure 1 
depicts such a case. We solved the one dimensional Schrodinger equation for a packet as 
described in |] , with and a barrier of very large strength. We assumed a total number of 
particles of N = 5 * 10 21 , obtained from the tabulated density of liquid helium and a volume 
of 1 cm 3 . We used a plate (potential) of width w=l cm and strength Vq = 4 eV ~ 10 15 , 
in units of % = 1. We took a detctor size of dx = 1 mm. It appears that, due to the large 
number of particles to be detected a simple weighing technique might be feasible, in this way 
we also avoid interaction of the packet with the detector. 

It is perhaps too optimistic to expect the proposed experiment will work so cleanly as in 
the simulation, due to all kinds of effects inside the drop that lead to decoherence as well as 
production of internal excitations like rotons, vortons, etc. However, some remainder of the 
effect might still show up in the counter. If it does, it will be a triumph for the quantum 
mechanical description of matter waves by means of wave packets of macroscopic size. 
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Figure 1: Number of Helium atoms counted at a distance of 5 cm behind a plate of thickness 
w=l cm as a function of time in seconds 
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